Casimir force in the rotor model with twisted boundary conditions 
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We investigate the three dimensional lattice XY model with nearest neighbor interaction. The 
vector order parameter of this system lies on the vertices of a cubic lattice, which is embedded in a 
system with a film geometry. The orientations of the vectors are fixed at the two opposite sides of the 
film. The angle between the vectors at the two boundaries is a where < a < n. We make use of the 
mean field approximation to study the mean length and orientation of the vector order parameter 
throughout the film — and the Casimir force it generates — as a function of the temperature T, the 
angle a, and the thickness L of the system. Among the results of that calculation are a Casimir 
force that depends in a continuous way on both the parameter a and the temperature and that can 
be attractive or repulsive. In particular, by varying a and/or T one controls both the sign and the 
magnitude of the Casimir force in a reversible way. Furthermore, for the case a = 7r, we discover an 
additional phase transition occurring only in the finite system associated with the variation of the 
orientations of the vectors. 

PACS numbers: 64.60.-i, 64.60.Fr, 75.40.-s 



I. INTRODUCTION 



For an 0(n),n > 1 model of a d-dimensional system 
with a temperature T and geometry oo d_1 x L the ther- 
modynamic Casimir force is defined by pj, [2[ 



p(r) 
Casimir 



(T,L) 



df£\T,L) 



dL 



where fex (T,L) is the excess free energy 

f£(T,L)=fV(T,L)-Lf b (T), 



(1.1) 



(1.2) 



and the superscript t denotes the boundary conditions. 
Here f( T >(T,L) is the full free energy per unit area of 
such a system subjected to the boundary conditions r 
and fb is the bulk free energy density. Accumulated evi- 
dence (H-[39j supports the conclusion that if the boundary 
conditions are identical — or sufficiently similar — at both 

(t) 

surfaces bounding the system, ^Casimir wm ^ e negative- 
In the case of a fluid confined between identical walls this 
implies that the net force between those walls due to the 
fluid will be attractive for large separations. On the other 
hand, if the fluid wets one of the walls while the other wall 
prefers the vapor phase, then the Casimir force is repul- 
sive. This implies that if the boundary conditions differ 
sufficiently the Casimir force can be expected to be posi- 
tive, or repulsive, for the entire range of thermodynamic 
parameters. In the intermediate case in which one of the 
surfaces of a system belonging to the Ising universality 
class has a weak preference for one of the phases of the 
fluid while the other one exhibits a strong preference, or 
for a given ratio of the surface fields and/or of surface en- 
hancements on both surfaces, it has been recently demon- 
strated [H, H3, that one can observe much richer be- 
havior, with the Casimir force changing its sign once, or 
even twice [36| . as the temperature is adjusted. In addi- 
tion, in [Z(| it has been shown via Monte Carlo simula- 
tions that in a system with a geometry L± x L?, subject 



to periodic boundary conditions both the magnitude and 
the sign of the Casimir force depend on the aspect ratio 
a r = L±/L». In this case general arguments have been 
advanced to suggest that at the bulk critical point the 
Casimir force vanishes for a r = 1 and becomes repulsive 
for a r > 1. These results are supported by exact calcu- 
lations for the two-dimensional Ising model. For further 
information regarding the results currently available on 
the critical Casimir effect the interested reader is referred 
to general reviews 0, Q as well as articles devoted to spe- 
cific aspects of the critical Casimir force [H, |41H43| . 

The critical Casimir effect and the corresponding 
Casimir force discussed above are due to spatial restric- 
tions imposed on the thermal average of the order pa- 
rameter and on its long ranged fluctuations in a system 
undergoing a second order phase transition. Based on 
analogy, some hints regarding the general behavior of the 
thermodynamic Casimir force can also be extracted from 
the information available on the quantum Casimir effect 
which is due to the spatial restrictions imposed on the 
possible fluctuations of the electromagnetic field. A series 
of reviews devoted to different aspect of this latter effect 
are available [45l - l55| . For the quantum Casimir effect it 
has been demonstrated that if the boundary conditions 
t are symmetric so that reflection positivity holds, the 



Casimir force ^c asimir is attractive [32j, |33|. According 
to the theory of Dzyaloshinskii, Lifshitz, and Pitaevskii 
[5(| in any system with a slab-like structure in which a 
material B separates two identical half-spaces A = C the 
force is attractive as well. When B is a vacuum, again 
according to [56| . this remains true even when the half- 
spaces A and C are not identical. This prediction, up 
to now, has been verified for all materials for which the 
Casimir force has been measured. Theoretical predic- 
tions exist, however, suggesting that in the latter case a 
repulsive Casimir force can be generated by special selec- 
tion of the material properties of A and C; see, e.g., [57j . 
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However, such a situation has not been experimentally re- 
alized. The omnipresence of attractive quantum Casimir 
electromagnetic force for objects in vacuum or air affects 
the work of micro and nano-machines [Hi], [HU, [HH], [H|| 
and might cause sticking of their working surfaces. The 
possibility of realizing and controlling a repulsive critical 
Casimir force might be one of the ways of overcoming the 
above-mentioned difficulties. 

In an attempt to shed additional light on the influence 
of differing boundary conditions on the critical Casimir 
force, we consider a film system with oo d_1 x L geometry 
consisting of local dynamical variables, say magnetic mo- 
ments, possessing 0(2) symmetry and constrained to lie 
in the x-y plane. The moments in one of the bound- 
ing surfaces are constrained to point in the same di- 
rection in that plane and to be oriented at an angle a 
with respect to the similarly aligned spins in the other 
bounding surface. Furthermore, in the case in which the 
moments have variable amplitudes, those amplitudes are 
fixed at a non-zero value. Alternatively, one might think 
of the studied system as a lattice gas of elongated, say, 
rod like molecules embedded on a lattice. We investi- 
gate both the equilibrium behavior of those moments (or 
molecules) — and the Casimir force that arises as a re- 
sult of that behavior — as a function of T and a. As 
we will see, among the results of our calculation are a 
Casimir force that depends in a continuous way on both 
the parameter a and the temperature and that can be at- 
tractive or repulsive. In particular, by varying a and/or 
T one controls both the sign and the magnitude of the 
Casimir force in a reversible way. 

We will refer to the boundary conditions described 
above as "twisted" boundary conditions. Subject to 
them, the moments within the system settle into a state 
in which they rotate with respect to each other as the 
region between the boundaries is traversed, creating a 
diffuse interface within it. The normalized excess free 
energy per unit area of the system, /J" , can be related 
to the corresponding quantity in a system with zero twist 
of the moments, which we will term a system with (+, +) 
boundary conditions. The quantified version of this rela- 
tionship is 

f£\T, L) = f£' + \T, L) + ^-a 2 TW(T, L), (1.3) 

where T^(T,L) is the finite -size helicity modulus [34 , 
IHH [6fl | that characterizes the energy of the system re- 
lated to the diffuse interface in it. The excess free energy 
/ix ^ (T, L) can also be resolved into regular and singular 
parts: 

/&>(T, L) = /W reg (T, L) + f£l sins (T, L). (1.4) 

In the case of the singular part of the excess free energy in 
the vicinity of the bulk critical point one has, according 
to finite-size scaling theory Q , 

fS si n g (T,L) = L-^X^( Xt ) (1.5) 



As a consequence, neglecting the "background" contri- 
bution to the Casimir force, one obtains 

F^(T,L) = L-*X<&(x t ). (1.6) 

Here Xt = a t tL x l v is the temperature scaling variable, 
t = (T — T c )/T c is the reduced temperature, at is a 
nonuniversal scaling factor, while X^ and Xq^ s are uni- 
versal (albeit geometry-dependent) scaling functions and 
v is the corresponding (universal) scaling exponent that 
characterizes the temperature divergence of the bulk two- 
point correlation length, £, when one approaches the bulk 
critical temperature from above, i.e. £(t — > + ) ~ £o~i _1/ 
with £q~ being some system dependent metric factor. For 
the behavior of T( Q )(T,L) near T c from (fOJ) and (Tj~5|) 
one derives 

T (a) ( T; L ) = £,-(«J-2)jf &») f ( L7 ) 

where X^ again is universal scaling function. Requir- 
ing a L-indepcndent behavior of T in the limit L — > 00, 
one obtains T(T) = lim L ^oo T^(T,L), with T(T) > 0, 
X { * ] (x t ) ~ \x t \ {d - 2)v and, thus, Tft) ~ \t\^ d - 2 > , which 
is in a complete agreement with [6fJ. 

When such a diffuse interface is present within the sys- 
tem and T < T c from Eq. (| 1 . 3[) it is easy to see that 

F&LtoF < Tc) * \a 2 T{T)L- 2 , L -+ 00. (1.8) 

Since T(T) > 0, Eq. (|1.8[) implies that the Casimir force 
will be repulsive and, for d > 2, much stronger, of the 
order of L~ 2 , than in systems with a compact interface 
where it is either of the order of L~ d , or smaller. 

The structure of the article is as follows. In the next 
section[TT]we define a lattice three-dimensional mean-field 
XY model and present numerical results for the behav- 
ior of the Casimir force within it. Section IIIII presents 
analytical results for the scaling function of the Casimir 
force within the Ginzburg-Landau mean-field theory of 
the three-dimensional XY model. In both sections HIl and 
IIIII we find interesting behavior of the force at a tempera- 
ture Tkink below the critical one of the bulk system when 
a approaches ir. We study this special case in section 
IIVI We deduce the existence of an additional second- 
order phase transition that is specific to this finite sys- 
tem. The article closes with a discussion presented in 
section [V] Technical details of the derivations are pre- 
sented in appendixes at the end of the article. 

II. THE CASIMIR FORCE IN THE LATTICE 
THREE-DIMENSIONAL MEAN-FIELD XY 
MODEL 

Consider a lattice of dimensions oo d_1 x L, with each 
site populated by an 0(2) fixed-length magnetic moment 
of magnitude m. We split up the lattice into (d — 1)- 
dimensional planes labeled 1, . . . , N, where L = Na, with 
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a being the lattice constant taken in the remainder to be 
equal to one. By translational invariance and neglecting 
the fluctuation within the planes, all moments in plane 
i must take the same value, equal to their mean value, 
and must point in the same direction. However, due to 
the anisotropy along the finite dimension, the moments 
will vary between planes. Let the moment in plane i be 
m.i. We take a nearest-neighbor coupling with strength 
J both in the plane and out of it, and so the energy of a 
moment in plane i will be 



Ui = — Jnii ■ (2(d — l)m.j + m, x + m i+1 ) 



(2.1) 



or, defining an effective magnetic field H; = J{2(d — 
l)nij + + m i+1 ) at that site, L7j = m t • H;. 

Approximating the moment as being isolated, in 
an external magnetic field Hi, we can assign it the local 
partition function 



Zi = 



iHi 



= 2irI (PmHi) 



(2.2) 



with j3 — (ksT) 1 and 9 the angle between Hi and m,-. 
On average, the component of rtii along Hi is 

ldbxZ hiPrnHj) 
(mcosfl) = — — = m^^r. — ttt, (2.3) 



P dH t 



where Iq and I\ are the corresponding modified Bessel 
functions of the first kind, while the component normal 
to Hi is 

1 r 2n 

(m sin 9) = - / d9msin9 e 0mHiCOs6 = (2.4) 
Z Jo 

so that the averaged moment is entirely along H;. Insert- 
ing the definition of Hi in terms of neighboring moments, 
we see that mi must satisfy the equation 



2(d - l)m.j + m. t _! + m i+ i 
\2(d - l)mj + m,-i + m,+i| 
x R (pmJ \2(d - l)mj + mi_i 



m 



i+i 



|)m (2.5) 



in the mean-field approximation (all quantities now im- 
plicitly averaged), with R(u) = h(u)/ Io( u )- 

If we define mo = m^r+i = for notational conve- 
nience, then the function 



N 



f({m l },N) = J2 



-m i -H i --ln (Io(pmHi)) 



(2.6) 



may be regarded as the total free energy functional of 
the system, because minimizing with respect to {m.,} 
yields the self-consistency conditions (|2.5p (See Appendix 
|A"]) . In order to compute the Casimir force on the sys- 
tem, we must also find the free energy per site in the 
bulk, i.e. when the system is very large. In that case, 
the moments {mi} will all be identical (at least near 
the center) and Eq. (|2.5[) tells us that the equation 



w^buik = Tifi(2dm/3Jmi j buik) determines their common 
magnitude m^buik- For T > dJm 2 /fee, the only solution 
is ^i,buik = while for T < dJm 2 /ks, there is a non- 
zero solution. Thus, the bulk system exhibits a transition 
at Tc.buik = dJm 2 /ks between an ordered (T < T c .buik) 
and a disordered phase. 

Letting a be the lattice spacing along the finite dimen- 
sion of the system, the bulk free energy density is 



dJm^huWL - -3 m {h (2dm/3Jm^buik)) 



and the Casimir force will be computed as 



(2.7) 



Fc 



mn — qj^ (/ L fb) 



A-Mittl. (2 . 



We now restrict our attention to the case d = 3 and 
investigate the system numerically. This amounts to solv- 
ing the simultaneous Eqs. (|2.5[) subject to particular 
boundary conditions on mj and mjy- We are interested 
in how the behavior of the system depends on the twist 
angle, a. 

When the twist angle is zero, the Casimir force is 
purely attractive (i.e, negative), as expected for match- 
ing boundary conditions. These features arc illustrated 
in the plot of Casimir force versus reduced temperature 
(Fig. 2). As the twist is increased, a low temperature 
region of repulsive Casimir force emerges. At a twist an- 
gle of 7r/2, the Casimir force becomes purely repulsive. 
When the system nears anti-symmetric boundary con- 
ditions, a = 7r, the Casimir force develops a kink at a 
temperature T kink < T c , bu ik- 

The nature of the kink will be discussed at greater 
length in a subsequent section, but we begin to under- 
stand it from the renderings (Fig. [T]). Below the kink 
temperature, the moments achieve the twist of tt by ro- 
tating about the axis while maintaining almost their full 
length m. Above the kink temperature, the moments all 
reside in a plane, and the twist is localized to the cen- 
ter of the system, where the magnetization has shrunk to 
zero. The nearest-neighbor interaction between moments 
imposes a free energy penalty for both rotating with re- 
spect to neighbors and varying in length. The transition 
indicates the point at which these penalties trade off in 
dominance. 



III. THE CASIMIR FORCE IN THE 
GINZBURG-LANDAU MEAN-FIELD THEORY 
OF THE THREE-DIMENSIONAL XY MODEL 

We now consider the continuous analogue of this sys- 
tem in the Ginzburg-Landau mean-field theory. The or- 
der parameter of the system is the magnetization profile 
m(z), where z is the finite dimension of the system. The 
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FIG. 1. (Color online) Renderings of the moments (d = 3, N = 30, a = n) for temperatures above (left) and below (right) the 
temperature where a kink occurs in the Casimir force. 



behavior of the system is found by minimizing the free 
energy functional (per unit area), 



,L/2 


~b 


dm 


/ dz 




l-L/2 


2 


dz 



-at Iml 2 
2 1 1 

+ 19 l m l 



, (3.1) 



with respect to m and subject to certain boundary con- 
ditions. The quantity t represents the reduced tempera- 
ture. 

Switching to polar coordinates, 



m(z) — ($(z) cos (p(z),$(z) sin <p(z)) 
the free energy functional is rewritten as 

r-L/2 
l-L/2 



(3.2) 



,L/2 

T[$,(p;t, L\ = dz 

J -LIT, 



dz 



I \ dz 
1 



-at®' 



1 



,g$ 4 



(3.3) 
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FIG. 2. (Color online) The Casimir force in the lattice model 
(d = 3, N = 30) as a function of reduced temperature t = 
(T — T Cjb uik)/r Ci buik for various values of the twist angle, a. 



Minimization with respect to tp(z) gives 
d 



which leads to 



-l*V]=o 



(3.4) 



(3.5) 



with an integration constant P v , independent of z, which 
roughly indicates the degree of twisting in the system. 
The condition from minimizing with respect to is sim- 
ilarly computed as 



d(f 
dz 



at$ + g§ 



or, with the identification (|3.5p . 

at<S> 



g<f> 3 . 



(3-6) 



(3.7) 



The problem is now that of solving Eq. (|3.7j) subject to 
twisted boundary conditions: 



¥>(±L/2) = ±a/2, 
$(±L/2) = oo, 



(3.8) 



i.e. where the moments at the boundaries are twisted by 
an angle a relative to one another. 

Note that, because of reflection symmetry in Eq. (|3.7[) 
and the boundary conditions imposed on <£>, we have that 
<J>(z) = z) and, thus, $'(z) = — $'(— z), whence 
$'(0) = 0. From the symmetry of Eq. (|3.5j) one con- 
cludes <^(z) = — (p(— z) which leads to ip(0) = 0. 

Multiplying Q3.7P by d$/dz and integrating with re- 
spect to z, we find a first integral 



P* = -\b 



P- 

_£ 

$2 



+ 



/d$ N 

\ dz 



1 



1 



of $ 2 + -g$ 4 , (3.9) 



with P$ being another integration constant independent 
of z. Let $o = ®{z = 0) be the amplitude of the order 
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parameter at the center of the interval. Then, taking into from which it follows that 
account that $'(0) = one can conveniently express P$ 

as /rf$\ 2 9 /ll 



dz ) v \&o $ 2 



g -(^-^t), (3-11) 



where 



a = f, ff = f. (3.i2) 



p2 i i The last result allows us to express the boundary condi- 

j 



P-s- = -7; bZ £ + 77 ai$ n + T9 $ n, (3-10) tions as 



r roc i poo 

dz= / d$^ = d<S> - = (3.13) 



l L/2 ^difi n f°° d$ 1 

P|($^ 2 -$-2) + ^(<f>2 -#§) + § ($4 _ $4) 



and 5=/ «b^=P,/ ^ ; , (3.14) 



These equations relate the integration constants, P v and 
$0i t° the system's external parameters, L and a. 

The stress tensor operator for a system with the free 
energy functional (|3.ip is [6l[ 



Tu.i = b 



dm dm 

dx k dxi 



d-2 



$' 2 + $y 2 " 

d 2 



2 4 y 



O 5 d 



dxhdxi 



- S k ,V 



_4(d- 1 

Calculating the (T z z ) component, one obtains 

2 D 2~ 



$ 2 . (3.15) 



1. 



dz 7 $ 2 



i«^ 2 -i ff $ 4 (3.16) 



which, according to the general theory, is a z-indcpendcnt 
quantity equal to the pressure —df^ T \T, L)ldL between 
the plates confining a fluctuating medium [3, |6l| . In our 
case, we see that 



-Pet 



(3.17) 



by p. 91) . From Eqs. ()3.10j) and (|3.16l) one observes that 
the Casimir force (the excess pressure over the bulk one) 
in this system is 



Casimir 



(t,L) 



l.if 



2 ^ 2 



(3.18) 



where 6{x) is the Heaviside step function. Here, we have 
taken into account that the bulk free energy density fb 
for the system is f b (t < 0) = -(at) 2 /Ag while f b (t > 
0) = 0. As expected, the Casimir force is sensitive to the 
boundary conditions, through the quantities P v and $o. 

It is easy to show that i 7 casimir(^ L) obeys the expected 
scaling. Indeed, in terms of the variables 



P v = -zX v L- A , at = x t L- 2 (3.19) 

the Casimir force reads (note that f, Xt, etc. are all 
dimensionless) 

^Casimir (t, L) = ^L^X^Xt), (3.20) 

where 



X 2 /X 2 -X 2 (x t +X 2 ), x t >0 
X 2 /X 2 -a Xt +X 2 ) 2 , x t <0 



(3.21) 

Taking into account that mean-field theories for short- 
ranges systems are effective d = 4 theories, one con- 
cludes that Eq. (|3.20[) is in full agreement with the ex- 
pected scaling behavior (|1.6[) of the Casimir force. From 
Eqs. f|3 . 20[) and (| 1 . 8[) one can derive the low temperature 
asymptotic behavior of the scaling function X(^(xt) of 
the Casimir force. We find 



1 



-co, 



(3.22) 



6 



with Xc a constant. Appendix [Cl contains a derivation 
of the expression in Q3.22[) . In that Appendix, we obtain 

the asymptotic expression for X^^(xt) 



400 



4 q 2m - 1* 2 



4^21 



3a z 



(3.23) 

when xt — > — oo. Note that Eq. (|3.23|l implies that 
Xq = 1 for the XY mean-field model. 

One can further simplify (|3.21[) by introducing the con- 
venient combinations of scaling variables 

r = x t /Xl and p = X v /Xl (3.24) 

Then the scaling function of the Casimir force reads 

T > 



y w w _i x V-(i + t)L 

Casl ' ~ 1 A V - (1 + r/2) 2 ], r < 



Eqs. P~T3"|) and (|3~T4]l then become 

da; 



A 



l ,J(x - l)[ir 2 + x(l +t)+ p 2 } ' 



and 



,1/2 

a = 2pA 3 / 
Jo 



xKCY 



(3.25) 



(3.26) 



(3.27) 
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FIG. 3. (Color online) A plot of the dimensionless scaling 
function for the Casimir force, Xcas, versus Xt (also dimen- 
sionless) for several values of a. The dotted curve is the 
Casimir scaling function in the Ising-like case of a critical 
fluid under (+, — ) boundary conditions. For Xt above a cer- 
tain value, this curve coincides with the one for the model 
studied here when twisted by an angle a w tt. 



In order to determine the Casimir force scaling func- 
tion X^"J s (xt), all one needs to know is the behavior of 
Xq = Xo(xt\a) and X v (xt\a) = pX$ as functions of 
Xt = tXq at a given fixed value of the angle a. For that, 
one has to solve Eqs. (|3.26[) and (|3.27l) after determining 




FIG. 4. (Color online) The scaling functions for the Casimir 
force in the Ginzburg-Landau model (solid curves), overlaid 
with those from the lattice model with N — 50 (data points) , 
for several values of a (from top to bottom: a — 0.987T, a = 
2tt/3, a = tt/2, a = tt/3, a = 0). 



the function from 



C-5 



dx 



2X o J[Xs,(C)/x ]i V(x - l)[x 2 + x(l + t) + p 2 } 

(3.28) 



with < C < 1/2, which directly follows from (|3~TTj) . The 
detailed knowledge of the behavior of the phase angle pro- 
file ip(() is not needed. The analytical treatment of Eqs. 
(|3.26[) - p.28[) is performed in Appendix iBl The numerical 
evaluation of the expressions derived there leads to the 
results for the Casimir force presented in Fig. [3] The 
comparison shows excellent agreement between the con- 
tinuum and the lattice model results from the previous 



1500 



1000 



8 



500 




FIG. 5. (Color online) Casimir force curves in the Ginzburg- 
Landau model for several values of a, overlayed with their re- 
spective asymptotic expressions (dotted) given by Eq. (|3.23|) 
as proven in Appendix C. 
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FIG. 6. (Color online) Left: The positions x[°q of the zeros of the Casimir force in the (x t , a)-plane. Right: The Casimir 

amplitude A^' s as a function of the twist angle a. This curve has been initially reported in [?| based on representation (|3.29[) 
derived there. The Casimir amplitude changes sign at a — n/3. 



section - see Fig. HI In order to demonstrate it, we scale 
the lattice results (t, F Cas imir) to (a t N 2 t, a F N 4 Fcasimir), 
where the scaling factors at and ap are determined by 
forcing the Casimir force with a = to agree between 
the two models. This was done numerically for N = 50, 
where we find a t w 2.977 and a F w 7.480 x 10" 5 . Fig. 
[5] shows a comparison between the low temperature be- 
havior of the Casimir force with the analytically derived 
asymptotic behavior reported in Eq. (|3.23p . We find 
that, for all cv, the asymptotic behavior is achieved for 
x t < -150. 

From Eq. p.25[) one can also infer some general prop- 
erties of the Casimir force. Taking into account that, at 
any fixed Xt and a, p is a definite function of xt and a, 
i.e. that p = p(r\a) one can, e.g., determine the coordi- 
nates x" of the zeros for the Casimir force for a given 

angle a. According to Eq. (|3.25[) one has that = 

for p(r\a) = \Jl + r, with r > and for p(r\a) = l + r/2 
when —2 < r < 0. A plot of the positions of these zeros 
in the (xt , a)-plane is presented in Fig. [5] 

The figure demonstrates how, by changing, e.g., the 
twist angle a, one can, at a given temperature t, make 
the Casimir force either repulsive or attractive. For 
< a < 7r/2, this can also be achieved by changing 
the temperature, i.e. the scaling variable xt, at a given 
fixed value of a. We also conclude that, when a — > 0, the 
position of the zero value of the Casimir force approaches 
— oo. This implies that when a = 0, the Casimir force will 
be attractive for all temperatures. Actually, for a = 0, 



X, 



Cas( x *) coincides with the known result for the Ising 

model system Q • The behavior of ^cts + ' ( x ) ^ s sri0wn as 
a thick black line in Fig. [3] These results are briefly re- 
derived in Appendix |B1 for the convenience of the reader. 

7r/2, we observe in Fig. [S]that x 



When a ~^ tt/2, we observe m t'ig. [OJtnat x' t ' — > oo. 
Thus a > tt/2 implies that the Casimir force will be 
repulsive for all values of Xt- As a increases, the repulsive 



force becomes stronger. We see from Fig. [3] that, when 
a = tt/2, the force is practically zero for all temperatures 
above the critical temperature of the finite system, while, 
for a > tt/2, it is repulsive in the whole temperature 
region. The cases a = 2ir/3 and a = 0.987T illustrate 
these features in the figure. One observes numerically 
that, for xt > —10, the a = 0.98-7T curve agrees with 
that of the mean-field Ising model with (+, — ) boundary 
conditions. At lower temperatures, there is an abrupt 
departure from the Ising model which will be discussed 
in the section to follow. The analytical expressions for 
the Ising model are known from [jj- For completeness, 
these results arc recalled in Eq. (|B23|) of Appendix [5] 



The behavior of the critical Casimir force, A 



(a) 
Cas 



XcJ s (x t = 0)/3, as a function of a, is illustrated in Fig. 
[SI Note that the Casimir amplitude becomes zero at 
a = 7r/3, so that the Casimir force for a = 7r/3 changes 
its sign at xt = 0. This was initially reported in Q and 
may also be seen in Fig. [HI Note that in Q different type 
of parametrization of the amplitude and phase profiles is 
used — they are parametrized via the Casimir amplitudes. 
This led to a restriction of the results presented to the 
critical temperature only. For example, for the determi- 
nation of the Casimir amplitudes in Q one has to solve, 
in our notations, the following system of equations (see 
Eq. (3.16) in 0) 



and 



[x 3 1 + (x l)X^A^ s ]-^dx (3.29a) 



i + av 4 a£> 



(3.29b) 



-V-l + (*-W 4 A&]- 1/a d* 
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IV. THE TRANSITION AT a = n 

The case a = tt warrants special investigation because 
it features behavior reminiscent of a phase transition. As 
mentioned in Section 3, the high temperature behavior 
of the system at a « tt tracks that of the Ising model. 
However, we find that a kink develops in all quantities 
in the system at a temperature tkink below the bulk crit- 
ical temperature of the system, and the system changes 
its character at this temperature. The lattice model also 
featured such a kink. In Section 2, we illustrated how the 
lattice system switches from a "rotational" state below 
the kink temperature to a "planar" state above it. We 
note that this phase transition-like behavior exists only 
in the finite system under the given boundary conditions 
and not in a thermodynamic sense. Below the transition 
temperature, the two coexisting phases are the rotational 
states with rotation plus or minus tt. There is sponta- 
neous symmetry breaking when the system orders in one 
of them. Above the critical temperature, there is a single 
state the "planar" one. 

The system incurs free energy penalties when adjacent 
moments vary in length or direction. We see how each 
type of state mentioned would extremize the energy. The 
moments in the planar state minimize rotation: they re- 
side in a plane and shorten to a length of zero at the cen- 
ter of the interval, where an abrupt reversal of direction 
occurs. The moments in the rotational state minimize 
length variation while gradually rotating from one end of 
the interval to the other. The temperature at which the 
kink occurs is the point at which the two energy penal- 
ties trade off in dominance. According to this descrip- 
tion, the planar state is characterized by Xq = X^ = 0. 
Indeed, this is what comes out of the Ginzburg-Landau 
model above the kink temperature (see Fig. [7|). While 
these quantities vanish at high temperature, their ratio 
Xip/Xo remains non-zero at all temperatures. 




FIG. 7. (Color online) The dimensionless quantities Xq 
(solid), X^/10 (long dashing; scale reduced for ease of plot- 
ting) and X v /Xo (short dashing) as functions of xt when 
a ~ 7r. The vertical dotted line indicates Xt = Xt.kink- 

In order to determine the kink temperature, we solve 



Eqs. (|3 . 13[) and (|3.14[) simultaneously in the vicinity of 
that transition point, i.e. with Xq = and X v /Xq finite 
but unknown. We find numerically (see Appendix [D]) 



^i,kink 



-28.1099. 



(4.1) 



Additionally, series expanding the conditions (|3.13p and 
p,14[) . we can find a first approximation to Xo(xt) in the 
vicinity of the kink. The result is that 



X (x t ) OC (Xt.kink - Xt) 



1/2 



(4.2) 



which agrees with standard mean field results for, e.g. 
the magnetization of a ferromagnet. 

We now present another, more transparent, analysis of 
this transition. We claim that, at high temperatures, the 
system's only energy extremum is the planar state, while, 
at low temperatures, the system has access to the planar 
state, as well as two energetically equivalent rotational 
states with a gradual turn of either tt or — tt. As previ- 
ously described, the system favors the rotational state at 
low temperatures and ignores the planar state. Such a 
situation is commonly found in Ginzburg-Fandau mod- 
els, where the free energy has terms quartic and quadratic 
in a variable or field of interest. Depending on the coeffi- 
cient of the quadratic term, there will be one fourth-order 
minimum at the origin or else a maximum at the origin 
and minima elsewhere. The planar state can be either an 
energy minimum or an energy maximum, while the ro- 
tational states are the symmetry-breaking minima that 
appear at sufficiently low temperature. By employing an 
approximation, we will show that this description fits the 
present system. 

It is useful to turn to the scaling variables (|3.19[) , in 
which the free energy functional, Eq. (|3.3[) . becomes 



T 



1/2 



<K 



dXi 



K 

Xq 



b 1 

9 L3 j -i/2 

+x t Xg 

while the amplitude equation (|3.7[) reads 



d 2 Xj 



—|- + x t A$ 



2Xi 



(4.3) 



(4.4) 



Due to the boundary conditions Q3.8P the amplitude pro- 
file near the edges, A$(£ ss ±1/2), is largely temperature 
independent, i.e. it is almost the same above and below 
the transition temperature. The behavior near £ = 
must therefore account for the physics of the transition. 
Near the transition temperature, we have Xq w and 
thus X*(C) < 1 for C ~ 0. Then Eq. ([44]) reduces to 



d 2 X<s, 
d( 2 



K 
X$ 



with solution 



**(C) = 




(4.5) 



(4.6) 
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denning M = X v /X . 

The free energy integral can now be computed in 
closed form using the asymptotic expression (|4.6j) . The 
expression for X<s>(Q is only valid up to some cutoff 
C = Y < 1/2, but the free energy from the edges of 
the interval (|C| > Y) will not contribute to the behav- 
ior of the system, provided that Y is large enough. The 
result of the integral is 



T 



2bY 
W 



xi 



-Y 
3 



-Y 4 M A 
5 



-x t Y 2 M 2 + M 2 
3 



(4.7) 



which is quartic in Xq. Recall that X$ is the (scaled) 
amplitude of the order parameter at the center of the 
interval. In the planar state, Xq = 0, while rotational 
states have X > 0. Therefore we expect to always see an 
extremum at X$ = 0, corresponding to the planar state, 
which will be a free energy minimum at high temperature 
and a maximum at low temperature. When the planar 
state is a maximum, two minima (rotational states) with 
\Xq\ > should emerge. Indeed, this behavior is clear 
from the form of (|4.7[) . The position of the non-zero 
minimum is found to be 



X = ±J-Ux t + ^M 2 



(4.8) 



which is only real for sufficiently low temperature: x t < 
-2Y 2 M 2 /3 < 0. The transition occurs when equality 
holds. This does not fully determine the temperature of 
the transition because both Y and M are functions of 
x t . The additional constraints arc afforded by matching 
the hyperbolic expression (|4.6|) for X$(Q with another 
expression correct near the edge of the interval. 

When C ~ ±1/2, A$ -> oo and, according to (|4.4j) the 
amplitude profile is determined by 



dC, 2 



XtX® + 2Xl 



which, with x t < 0, is solved by 



**(0 



(4.9) 



(4.10) 



Now imposing the continuity of X$(£) and X'^(Q from 
(14.6[) and (|4.10l) at £ = Y allows us to solve for the pa- 
rameters at the transition point. Proceeding numerically, 
we find: 



as (x^kink — Xt) 1 / 2 , in agreement with the power law pre- 
viously found. Finally, it is easy to show that at the kink 
temperature the rate of the change of the phase in the 
middle of the system diverges. Indeed, from Eq. (|3.5[) 
and using the definitions (|3.19[) one obtains 



x*(0 



Xn 



MO _ x. 



X 2 



M 
X 



= = — . (4.12) 



Thus, in the limit £ — > one derives at xt = Xt,kink 
that (d(p(()/d()\^ =0 = M/Xq — > oo, since A 0i kink = 0. 
Therefore, at the kink temperature all the change in the 
phase of the moments happens at the middle of the sys- 
tem, where there length becomes zero. The phase of the 
moments jumps there from to tt. 



V. DISCUSSION AND CONCLUDING 
REMARKS 

We have studied the 0(2) model in a three-dimensional 
film geometry, and find identical predictions from lat- 
tice and continuum mean-field theories. Consistent with 
systems of similar type, the Casimir force with symmet- 
ric boundary conditions (a = 0) is attractive, while the 
Casimir force with anti-symmetric boundary conditions 
(a = tt) is repulsive. In particular, the critical Casimir 
force, i.e. -Fcasimh- at the bulk transition temperature, 
changes from repulsive to attractive. This is a stan- 
dard result, but we also find intermediate scenarios when 
< a < tt which feature critical Casimir forces, and scal- 
ing functions for the Casimir force (illustrated in Fig. [2]), 
different from those of the symmetric and anti-symmetric 
cases. The Casimir force may therefore be continuously 
adjusted at constant temperature by varying the twist a, 
or at constant twist by varying the temperature. 

Additionally we find that, when the boundary condi- 
tions are perfectly anti-symmetric (a = tt), the system 
undergoes a phase transition at a temperature below the 
bulk critical temperature. We are able to understand 
this transition as a symmetry-breaking effect: at high 
temperatures the moments of the system are confined to 
a plane, while at low temperatures they rotate about the 
z-axis by either tt or — tt to satisfy the boundary condi- 
tions. The high temperature behavior tracks that of an 
Ising model whose order parameter is always simply up 
or down, but our system departs from that behavior at 
the transition point, once moments find it energetically 
favorable to rotate. 



zt,kink w -22.4587, M kink ~ 19.4498, 

*o,kink = 0, r kink PS 0.2984. (4.11) 

Compared to the exact numerical results obtained in Ap- 
pendix [D] these values are consistent as a first approx- 
imation, as are the behaviors of Xo(xt) and M(xt). In 
particular, the leading order contribution to Xo(xt) goes 



ACKNOWLEDGMENTS 

D. Dantchev acknowledges the partial financial sup- 
port of Bulgarian NSF, grant No. DO 171/08. J. Rud- 
nick acknowledges partial support from the NSF through 
grant DMR-1006128. 



10 



Appendix A: Lattice free energy 



We claim that 



N 



f({ mi },N) = J2 



im 4 -Hi-i In (Jo (pmHi)) 



(Al) 



is the free energy functional of the lattice model con- 
sidered in Section 2. Indeed, here we will demonstrate 
that minimizing it with respect to the leads to the 
mean-field consistency equations 



Appendix B: The amplitude and phase profiles and 
Casimir force derivation within the 
Ginzburg-Landau mean-field theory of the 
three-dimensional XY model 



In this appendix, we derive some analytical expressions 
needed for the numerical evaluation of Eq. (|3.2ip for the 
scaling function of the Casimir force. 



m, 



l HiI {fimHi)' 



(A2) 



Recall that H; = J(2(d — l)rrij + rrii_i + mj+i) and that 
we take mo = mjv + i = for notational convenience. 
Differentiating, 



= V m J = 2J(d-l) 
+ J mj_! — ? 
+ J 



HihiPmHi) 



Hi Io(PmHi) 





X-,X+<1 


r=-l+2p 







x-, x+e C 


-5 








X-<1,X+>1 1 


^JtT _ • 1 


-10 







Hi+i h(pmH i+1 ) 



(A3) 



H l+ i I (f3mH i+ i 

for i = 2, . . . , N — 1. If i = 1 (i = N), we find the same 
condition with the second (third) term omitted. Defining 

g, = m, — m -, 

HiloiPmHi)' 

we may write Eq. (|A3|) as 
= V mi / = 2J(d-l) gi 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
P 

FIG. 8. (Color online) The loci of points in the (p,r)-planc 
for which the roots x± possess the properties discussed in the 
main text. When any of the roots approaches the thick blue 
line, Xq —¥ oo. One observes that this is only possible for 



go = gw+i = 0, (A4) T < °> i- e - when rX o =xt^ -oo. 



Jgj_i + Jg 



i+i 



(A5) 



for each i = 1, . . . , N. In order to show that (|2.5[) holds, 
we must show that (|A5|) is only solved when gj = for 
all i. This is seen by writing the linear equations (|A5|) in 
matrix form, i.e. Ag = with tridiagonal N x N matrix 



We start by determining the behavior of the amplitude 
profile X&(Q - see Eq. ()3.28j) . In addition, wc will also 
determine the phase angle profile <p((). Note that, in 
terms of the scaling variables (|3.19p and (|3.24p , we obtain 
the phase angle ip(() as 



A = J 



( 2{d-l) 1 
1 2(d-l) 1 
1 







\ 



v(C) 



i 



(Bl) 



via Eq. (|33|) . Let 



(A6) 



x± 



-(r + l)± V(r + l) 2 -# 2 



(B2) 



V ••• 1 2(d-l)/ 

whose determinant is computed as 

J N 



dct A 



2^(d-l) 2 -l 

x 



N+l 



(d-i + ^{d-iy-i) 
-(d-i- v/(rf-i) 2 -i) 



JV+l 



(A7) 



which is non-zero for d > 2 and > 1, so the only 
solution is g^ = 0, as desired. 



be the roots of the quadratic term in the square brackets 
in the denominator of (|3.26l) . 

In order to perform the integration in ()3.26j) . where 
the integrand is a positive function for all points from 
the integration interval, one needs to know if these roots 
are real or complex (see Fig. 8). Thus, there are two 
subcases: it A) the roots are real, and B) the roots are 
complex conjugates of each other. 

First consider the subcase 

A ) The roots x± are real. 

In that case the positivity of the integrand implies 
X- < x+ < 1. Taking the above into account and us- 
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ing the corresponding expression reported in [62J , 
dx 



y V( x ~ - x+)(x - x-) 
-.F arcsin , 



1 -x. 



s/T=x^ \ 



y-x- 



x + — X- 

1-X- 



, (B3) 



provided y > 1 > x + > x_. From Eq. (|3.28[) . one 
obtains 



1 



1 



2 X a ^T^T 

x F 



1-x- 



arcsm , 



x + — X- 



, (B4) 



(Xj/Xq) 2 - a;_ ' V 
so that, for £ = with = 0) = Xq, it follows that 

2 



— X- 



1 



(B5) 



J 



and 



1 



2 



In (|B3[) and (|B4|). A refers to the elliptic integral of the 
first kind, while K in (|B5|I is the complete elliptic integral 
of the first kind. Solving (|B4| for Xj> one finds 



X|, 



c + 



1 



, (B6) 



where sn denotes the corresponding sine-amplitude Ja 
cobi elliptic function. Finally, inserting (|B6[) into (|B1| 
and (|3.27p and performing the integration, we arrive at 



x- - 1 
z_ - 1 



, am 



, am 



-x «yi - x-, 







x+ 


— X- 


1 




V i 


- x_ 






Vt 


— X- 








— X- 



, (B7) 
,(B8) 



where IL(n,<p, m) is the incomplete elliptic integral of the 
third kind and am(u, m) is the amplitude for Jacobi el- 
liptic functions. 

The relationship between x±, r and p is straightfor- 
ward. From (IB2I) we have 



- 1 — X— — x+ , 



P 



y/\x-x+\. (B9) 



Taking into account Eq. (|B5|) one can further simplify 
Eq. dEHJ to 



= s/\x-X+\X Q f 



=n 



i' 



1 - x_ 



(BIO) 



Here Xo, x_ and x + arc known functions of r and p. 
Fixing, for instance, r, one can solve (|B10[) numerically 
for p. Then, knowing Xq-, the scaling function for the 
Casimir force is found from Eq. (|3.25p . These scaling 
functions are plotted in Fig. [3] 
B) The roots x± are complex. 

In this case, the roots are complex conjugates of each 
other, i.e. X- = 3Tj7. Taking this into account and using 
and using the corresponding expression reported in [63j |. 



we obtain 



dx 



y tJ(x - l)(x - x+)(x - x-) 



= 7Tr F 



y — 1 — r 

arccos | | , w 

y — 1 + r 



where 



*(x-,x+) = - x-)(l - x+) 



T T" p 



and 



w 2 = w 2 (x-,x + ) = ^ + 



2 2y/(l-X-)(l-X+) 



1/ 1 3 + r 



2 V 2 yj2 + r + p 2 J 
According to Eq. (|3.28[) . the above implies that 



1 



2X ^¥ 



(X„/X Q f - 1 - r 
1 (X*/X )*-l + r l 



(Bll) 



(B12) 



(B13) 



(B14) 
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f f 



FIG. 9. (Color online) Plots of the amplitude profile and the angle of the order parameter for a = tt/3 and some choices of xt. 
We observe that, when the temperature increases, the value of the amplitude in the middle of the system decreases. The twist 
of the local variables through the system spans over the total system almost uniformly for low temperatures, while for higher 
ones it concentrates more and more in the middle of the system where the amplitude is at its smallest values. 



and, for £ = with X$(( = 0) = Xq, it follows that 



X = -=K (w) 
\ r 



(B15) 



Solving (|B14I) for X$ gives 



*l [ C + =: 



xa 



l-r- 



2r 



1 -cD.(2X oy /f(,w) 



(B16) 



where cn denotes the corresponding cosine-amplitude Ja- 
cobi elliptic function. Finally, inserting (|B16|) in (|B1[) and 
(|3.27[) and performing the integration, one arrives at 



r 



x <^ n 



l v 1 - w 2 
r — 1 



1 



\fuP- 



-n 



r - 1 
r + 1 



7T TV 

, - - am (2^X0 C, w) , , 5 = 
^ yr — 1 



pX 1 

C arccot 

l-r 2 



2Vfdn (2^o C»«0 



p sn(2VfX C, 



(B17) 



and, setting £ = in the above equation, we have 
pX 



a 



l-r 



2p_ 
2 - 1 



n 



r 

1 — w 2 

r-T 
r + 1 



n 



7T 

2 



(Vr~X , 



w 



'Si 



1 



(B18) 

where we have used that, according to Eq. (|B15I) . 
XoV^ 7 = 2K(u>), and that dn[2if(w),w] = 1 and 
sn[2K(ui), oj] — 0. The amplitude profile X$(() and the 
angle profile <p(£) are plotted in Fig. 9. Now, using the 
properties of the am and n functions, the above equation 
can be further simplified to 



a 



pX Q 



4p 



IV 1-w 2 



x n 



1 



w 



Vw 2 



1 



(B19) 



Recall that the relation of X- and x+ to r and p is given 
by Eq. (|B9[) . As in the previous subcase, X 0l x_ and 
x + are known functions of r and p, and this equation 
is solved numerically to produce the scaling function for 
the Casimir force. 

From the expressions derived above, it is easy to re- 
produce the results previously known for a = 0. As we 
will see, this provides a new representation of the older 
results which is quite convenient for numerical evalua- 
tion. First, let us note that, from Eqs. (|B8|) , (|B9[) and 
(|B18p . one immediately obtains p = 0. Thus, from Eq. 
(|B2p . it follows that we are in the subcase A) of real 



+ 



roots. Then x 

x + = — (r + l),x- 
find 



0,x_ = — (r + 1) for t > — 1 and 
= for t < — 1. From Eq. (|B5|) . we 



*o(t) 



2Jf U/-(r + i; 



/t + 2 



K 



T + l \ 
T + 2 J 



T < -1 
T > -1 



(B20) 



For the scaling function of the Casimir force, Eq. (|3.25[) 
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X 



(+,+; 

Cas 



(r) 



-4(r- 



-4 



2) 2 A 4 



AK 4 



(y/-(r + lj) , t < -1 

> T > -1 



T + l 

r + 2 



r+1 
(r+2) = 



A" 



(n/*D 



(B21) 



r > 



where we have denoted the a = boundary conditions as 
(+, +). Denoting the argument of the elliptic A function 

I 



I 

in a standard way with k and recalling that x t = tA 2 , 
the above expressions can be rewritten in the form 



-4 (i - /c 2 r a 4 (fc) , xt-- 

-4A 4 (fc), x t 
-4fc 2 (l - fc 2 )A 4 (fc) , x t 



-4{l + k 2 )K 2 {k), x t <-n- 
4(2fc 2 - l)A 2 (fc), 0>x t >- 
4(2fc 2 - l)A 2 (fc), x t >0 



(B22) 



The result (|B22|) was originally reported in [3]. The be- 
havior of g {xt) is shown as a thick black line in Fig. 
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r 



The scaling function of the Casimir force under (+, 
boundary condition in the Ising mean-field model is 



X, 



(+,-), 

Cas 



64fc 2 (l-fc r 
[2K(k) 
[2K(k)} 4 (l 



^ [K(k)] 

4 



2\2 



k 2 ) 



X t = 
Xt = 

x t 



-2 [2K{k)f {2k 2 
-2 [2K{k)f {2k 2 
2 [2K{k)f {k 2 - 



1), 
1), 



x t < 
< x t < 2tt 2 
x t > 2tt 2 



(B23) 



xj^l \x) is plotted (marked with filled circles) in Fig. 

m 



Finally, note that the scaling functions ^cas + ' ( x ) anc ^ 
(+,- 

Cas 



-X"cta J ust derived, are related through 



X { c t +) {x) = --X^'{-x/2), 



1 



-(+,-), 



(B24) 



and thus, for the corresponding mean-field Casimir am- 
plitudes, one has 



,(+,+) 

*Cas 
*Cas 



(B25) 



Appendix C: Derivation of the low-temperature 
asymptotic behavior of the Casimir force within XY 
Ginzburg-Landau mean-held model under twisted 
boundary conditions 



According to Eq. (|3.25|) when x t < 



X&>(r)=X V-(l + r/2) 



(CI) 



where r and p are defined in Eq. p. 241) . We need to find 

the behavior of X^^{t) for xt — oo. 

Let us first clarify what is meant by the asymptotic 
behavior of r and p in the regime x t — > — oo. For low 



I 

temperatures one expects $(z) — $>(z = 0) = $o and 
dip/dz ~ a/ A. From Eq. (|3.5[) and the definition given 



in Eq. (|3.19[) one then obtains X v ~ aX^ and, thus, from 
Eq. (|3.24[) . p ~ a/X(j. In terms of Xt the equation for the 
order parameter amplitude is given in (I4.4[) . Under the 
assumptions already made, the above equation becomes 
~ a 2 A + x t X + 2X$. One concludes that X > 1 
with 



An 



-{x t + a 2 )/2 



(C2) 



when x t — > — oo, and that p 2 + t + 2 ~ 0, i.e., that r — > 
— 2— p 2 when xt — > — oo. Thus, the regime which we need 
to consider in (jCip is r — > — 2 - 
For the Casimir force from Eq. 



9 2 

EH) 



with p ~ a/Xo <C 1. 



wc then obtain 



A^(* t ->-oo)^A V-p 4 /4]. 



(C3) 



It is easy to check that in the asymptotic regime of 
interest x± in Eq. (|B2[) are real. Thus, we need to study 
the asymptotic behavior of X given by Eq. (|B5[) , tak- 
ing into account the right-hand side of Eq. (|B10[) which 
relates Xq to a. 

Setting 



where p -C 1 and a 
becomes 



r = -2-p 2 + a, (C4) 
— ^ it is easy to show that Eq. 



X ~ ln[16/a], 



(C5) 



14 



while Eq. (jBlOl) simplifies to 
a = p(X 



(C6) 



2), 



where we have used Eq. (|C5|) . Note that Eqs. (|C5|) and 
(|C2[) imply that a is exponentially small in \/[x7[ and, 
thus, in the remainder we will omit a in Eq. (|C4|) and 
in any expansion that involves r. Expressing p from Eq. 
(|C6|) in terms of A and a and inserting the result in Eq. 
(|C3[) . one obtains an expression for the Casimir force in 
terms of Xq and a: 



x£ B (x t -> -co) 



a 2 (A 2 



4An 



12) . 

; 4 



(C7) 



Then, making use of Eq. (|C2[) , one obtains 



at that point. If we make the substitution A$ = Xqu, 
we see that f(M k , 0) is actually non-singular: 



f(M k ,0) = M k 



dy 



yVy 2 - 1 
l 



(A/ 2 + art.krt 2 + y 2 (y 2 + 1) A 4 ) 1/2 , 

G?y 7T 

= 2" 



(D4) 



This means that (df/dXo)(Mk,0) must vanish due to 
Eq. (|D3[) . The derivative is taken most easily from the 
expression in (|D4[) . giving 



i 



■4^/2^ 



3a 2 



(C8) 

where xt — > — oo. This is the result reported in Eq. (|3.23j) 
in the main text. 



Appendix D: Determining the kink temperature in 
the Ginzburg-Landau model 



We employ the scaled variable defined by Eq. (|3.19l) . 
When the boundary conditions arc anti-symmetric, we 
find a kink in the Casimir force at a temperature Xt ;k i nk 
below the bulk critical temperature. At that point, the 
two integration constants Xq and X v both switch be- 
tween being identically zero (xt > x*,kink) and being pos- 
itive (x t < Xt,kink)- Note, however, that the quotient 
M = X^/Xq remains non-zero for all temperatures. We 
determine the kink temperature by enforcing the bound- 
ary conditions, Eqs. (|3.13p and (|3.14|) . 

At the transition point, things are simplified because 
Xq, X v — > 0. The length condition (|3.13j) takes the form 



dXi 



{Ml 



kink 



X^kinkX 2 +A|) V2 ' 



(Dl) 



The twist condition must be treated with more care be- 
cause the integrand appears to be singular when Xq = 0. 
Without taking A to zero, (|3.14[) may be re-expressed 
as 



I- MX, 



dX, 



<1> 



X a Xq, yj A J, — Xq 

1 



(Af 2 +x 4 A 2 +A 2 (A 2 +A 2 )) 1/2 

= f(M,X ), (D2) 

which holds at all temperatures. In particular, just be- 
low the kink temperature, Xq is small but non-zero, and 
(abbreviating M k = M kink and x t . k = Zt.kink) 

| = f(M k , 0) + ^(M k , 0) X + 0(A 2 ) (D3) 



^-(M fc ,0) = -M k X 

O^Q Jl 



ydy 



x t , k + 2y 2 A 2 + 2A 2 



vV^T 

2 



(M 2 +x t , fc y 2 A 2 + y 2 (y 2 + l)A 4 ) 



4^3/2 



(D5) 



X o = 



Despite its appearance, this does not trivially vanish 
when Xq — > 0. Instead, restore the original variable 
A$ = X y to find 

9/ . , , f°° Acj> dX$ 



x t , k + 2A 2 + 2A 2 



(Af 2 + x t , fe A 2 + A 2 (A 2 + A 2 )) 



2\\3/2 



(D6) 



X =0 



which suffers no singularity when Ao is replaced by zero. 
Thus the second condition on xt k and M k is 



dA$ 

(A/ 



x t ,h + 2A 2 



x t , k X% + X 



4\3/2' 
<S>) 



(D7) 



Eqs. (|D1|) and (|D7[) may be recast, with the aid of Eqs. 
(|BT5|) and ([B19]) . into 



'M k = 2K 



(D8) 



and 
1 



M k y/x t ,k + 2M fc = AM k K 



(xtj. + 2M k )E 



X t ,k 


- 2M k 


X t ,k 


+ 2M k 


Xt,k 


- 2M k 


Xt.k 


+ 2M k 



, (D9) 



which are easily solved numerically to give 

x m nk « -28.1099, A/ kink « 21.5491. (D10) 

In Eq. (|D9|) . E(x) is the complete elliptic integral of the 
second kind. 
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